Semiclassical spin transport in spin-orbit coupled bands 
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Motivated by recent interest in novel spintronics effects, we develop a semiclassical theory of spin 
transport that is valid for spin-orbit coupled bands. Aside from the obvious convective term in which 
the average spin is transported at the wavepacket group velocity, the spin current has additional 
contributions from the wavepacket's spin and torque dipole moments. Electric field corrections to the 
group velocity and carrier spin contribute to the convective term. Summing all terms we obtain an 
expression for the intrinsic spin-Hall conductivity of a hole-doped semiconductor, which agrees with 
the Kubo formula prediction for the same quantity. We discuss the calculation of spin accumulation, 
which illustrates the importance of the torque dipole near the boundary of the system. 

PACS numbers: 72.10.-d, 72.15.Gd, 73.50.Jt 



Electrical control of spins in systems with spin-orbit 
interactions is of basic interest and has great poten- 
tial in semiconductor spintronics P, 0- Iri recent 
years, steady progress has been made towards realiza- 
tion of convenient semiconducting fcrromagnets and spin 
injection into semiconductors from ferromagnetic metals 
|j[ EJ S El The spin transport theory presented 
in this Letter is motivated generally by current interest 
in novel spin-related transport effects in semiconductors, 
and particularly by interest in various schemes that gen- 
erate spin-polarized currents m m e m m m un. 

Using a semiclassical wavepacket approach, we find that 
the spin current can be expressed as the sum of several 
physically transparent terms which are grouped together 
in a Kubo formula description. As an example, we use 
our theory to derive an expression for the intrinsic spin- 
Hall conductivity 0, |n| of a hole-doped semiconductor. 

Semiclassical formulations of transport theory exploit 
the smooth variation of transport fields on atomic length 
scales. Previous semiclassical theories of spin trans- 
port 0, 0, 0, 0, have not accounted explicitly 
for intrinsic spin-orbit interaction in the crystal apart 
from, occasionally, its role in the relaxation of non- 
equilibrium spin-polarizations. In this Letter, we ap- 
ply the wavepacket approach introduced by Sundaram 
and Niu |2l|, which captures the consequences of the 
wavevector dependence of the Bloch spinors, to treat spin 
transport in spin-orbit coupled bands. This wavepacket 
approach has already been successful in describing the 
anomalous Hall effect in ferromagnetic semiconductors 
[22| and transition metals (2^, interpreting it as a con- 
sequence of the Berry-phase correction to the group ve- 
locity induced by the intrinsic spin-orbit interaction. We 
show here that the Hall spin current in response to an 
electric field is non-zero even in paramagnetic systems 
and that, in addition to the Berry phase term evaluated 



in a recent paper , other contributions must be taken 
into account. First of all, there is a contribution from the 
electric field correction to the average spin orientation 
of a wavepacket. In addition, there are also contribu- 
tions from the spin dipole and torque dipole of a carrier, 
which arise from the fact that spin and torque distribu- 
tion within a wavepacket generally differ from that of 
the charge. Including all these contributions, we obtain 
a total semiclassical spin current which is in agreement 
with the Kubo formula expression for the same quantity. 
We show that nonequilibrium spin polarization near the 
sample edge is driven not by the spin current alone but 
by the sum of the spin current and torque dipole density. 

The semiclassical dynamics of each spin-charge car- 
rier in a non-degenerate band is described by a 
wavepacket, whose charge centroid has coordinates 
(r c ,k c ). Wavepacket construction is thoroughly ex- 
plained in [23, |24| . When expanded in the basis of Bloch 
eigenstates, the wavepacket has the form: 



(1) 



= / d A k a(k,i)e lkr |u(r c ,k,t)). 



In the above, the wavefunctions \u) contain [25j correc- 
tions linear in the electric field. They form a complete set 
and retain the Bloch periodicity. The function a(k, t) is 
a narrow distribution sharply peaked at k c , and its phase 
specifies the center of charge position r c . The size of the 
wavepacket in fc-space must be considerably smaller than 
that of the Brillouin zone. In real space, this implies that 
the wavepacket must stretch over many unit cells. 

In the presence of a uniform electric field, the semiclas- 
sical equations of motion for a non-degenerate band read 
El: 



Mr 



de 



eE x ft and Kk c = eE, 



(2) 



2 



where e is the carrier charge, s is the band dispersion, 
and ft is the Berry curvature of the Bloch state [2lj| . 
Henceforth, k c will be abbreviated to k. The effect of 
the electric field is thus twofold: it drives the center of 
the wavepacket in fc-space, and it gives rise to a non- 
adiabatic correction to the wavefunctions, which mixes 
the bands at each k. 

Following the strategy of Boltzmann transport theory, 
we consider a collection of particles described by a phase 
space distribution /(r c ,k, t). This distribution can drift 
according to the semiclassical equations of motion (2), 
and can also change due to collisions: 



df 

at 



df_ 



df 

dk 



4f, 
dt 



coll- 



(3) 



The collision term on the right hand side may be mod- 
eled by relaxation times or more accurately by collision 
integrals as usual. 

The spin density distribution is defined as 



S(r,t) = 



d 3 kd 3 r c f(r c ,k,t)(S(r-r)s), (4) 



where s is an arbitrary component of the spin operator, 
and the bracket indicates quantum mechanical average 
over the wavepacket with charge centroid (r c , k). Further 
analysis of this distribution will be facilitated by making 
the analogy with the standard coarse graining of electro- 
dynamics in material media |2fij| . Our wavepackets play 
the role of 'molecules', whose size will be taken as much 
smaller than the length scale of the distribution func- 
tion. We are thus allowed to view the 5-function in the 
above definition of the spin density as a sampling func- 
tion with a width somewhere between the microscopic 
scale of the wavepackets and the macroscopic scale of 
the distribution function. We can therefore write it as 
<5[(r — r c ) — (f — r c )] and expand it around r c , keeping 
only the first order term. Performing the integral over 
r c , the spin density can be re-expressed in the following 
form: 



S = / d a kf(s) 



d 3 kfp s 



(5) 



where / = /(r, k, t), and p s = ((? — r c )s) | rc=r is the spin 
dipole. The two terms can be regarded as monopole and 
dipole contributions. The second term is analogous to 
the contribution to the charge density in electrodynamics 
from the divergence of the polarization. 

Spin is in general not conserved, and for what follows 
it will be useful to define a quantity, which we shall call 
the torque density, in order to include the rate of change 
of spin into our discussion of transport: 



T(r,t) 



d 3 kd 3 r c f(r c ,k,t)(S(r-r)r), (6) 



in which f is understood as i[H, s] and H is the Hamil- 
tonian. Following the steps outlined above, the torque 



density becomes: 

T = J d 3 kf(f) - V J d 3 kfp 



(7) 



with the torque dipole p T = ((? — r c )f)| rc=r . 

We evaluate the spin-current using the microscopic 
spin-current operator and the semiclassical distribution 
function: 



3 s (r,t) 



d 3 kd 3 r c f(r c ,k,t)(5(r-r)is). (8) 



Throughout this paper, symmetrization of products of 
noncommuting operators is implied. After expanding, 
the spin current takes the form: 



J s = / d 3 kf(rs) - V • / d 3 kf((r - r)fs 



(9) 



For a homogeneous system, where the distribution func- 
tion is independent of position, the gradient term van- 
ishes, and it is permissible to use Bloch states (which 
may be regarded as the limit of very wide wavepackets) 
to evaluate the carrier spin current (rs). Since the Bloch 
states contain first order correction in the field, this can 
in general yield an overall linear-in-field spin current even 
with the equilibrium distribution function. This intrinsic 
spin current has been evaluated for a number of systems 
recently, and identical results are obtained with the semi- 
classical approach developed here. 

To illuminate the underlying physics, we decompose 
the carrier spin current into a number of terms: 



d 3 kf(r c (s) + 



dp s 
~dt 



Pi- 



(10) 



The first contribution is convective, arising from the 
fact that the total spin is transported as the wavepacket 
moves. The second comes from the rate of change of the 
spin dipole, while the third is from the torque dipole. 
This decomposition makes it possible to compare the 
Kubo formula result with those based on various heuris- 
tic arguments. The authors of restricted their scope 
to the convective term and considered only the Berry 
phase contribution to the charge center velocity r c . The 
present semiclassical decomposition allows us to recog- 
nize the missing terms due to the spin and torque dipoles, 
as well as a field correction to the carrier spin in the con- 
vective term. The approach of |ld| would give a zero 
result for the Rashba model, whereas the Kubo formula 
approach of which agrees with (TDJ, yields a nonzero 
spin Hall current for this model. Interestingly, the spin 
Hall current in the Rashba model can be obtained ex- 
actly from a heuristic argument based on a velocity and 
field dependent correction to the carrier spin as discussed 
in [TlT ]. This approach is, however, applicable only to 
single-band models with wavevector-dependent Zeeman 
coupling. 
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d z k%(s). 



(11) 



The spin density and current satisfy the following equa- 
tion of continuity: 

It is seen that the torque density appears in the source, 
accounting for the spin non-conserving terms in the 
Hamiltonian, and acting as a bulk mechanism for spin 
generation. The second term accounts for the effect of 
collisions. 

The source can be decomposed into intrinsic and ex- 
trinsic contributions, depending on the equilibrium and 
non-equilibrium parts of the distribution respectively. If 
we restrict our attention to homogeneous sytems the 
torque density is simply /(f). We find that this term 
is always first order in the electric field, and is given by 
/|E • ^j^-. We are thus justified in replacing / by its 
equilibrium value /o, in which case this term is purely 
intrinsic. The second term in the source, which depends 
on the nonequilibrium shift in the distribution function, 
is entirely extrinsic. 

Our formalism so far applies to independent non- 
degenerate bands, and for the Rashba model its predic- 
tions are in agreement with ^lj. There exists a parallel 
formalism for coupled degenerate bands, which yields the 
same results as given above )25| . In this case, the dis- 
tribution function becomes a density matrix, while (s), 
(f), r c , p s and p r are replaced by matrices. To find the 
macroscopic expectation values, one traces over the den- 
sity matrix. This formalism can be applied, for example, 
to the spherical four-band Luttinger Hamiltonian, 



Ho 



2m 



[(7i + 7T72)fc 2 - 2 72 (k-J) 2 ], 



(12) 



where J is the operator for angular momentum 3/2 and 
7i, 72 are the Luttinger parameters. The Bloch states 
are eigenstates of the angular momentum projection in 
the k direction, Jk- The four bands are split (for finite 
k) into two degenerate manifolds with Jk = ±3/2 (heavy 
holes) and Jk = ±1/2 (light holes). 

Let us take a closer look at the source term, using 
the four-band model as an illustration. This discus- 
sion applies to either of the heavy and light-hole man- 
ifolds. In equilibrium, the density matrix is diagonal 
and equal to fo for each band. The mean spin in the 
z-direction is (s z )± 3 / 2 — ±^f for the heavy holes, and 
it is (s z )±i/2 = ±insf for the light holes. The spin ex- 
pectation values have opposite signs in the two bands, so 
that, when averaged over the equilibrium density matrix, 
the intrinsic term in the source, /(f), will vanish. The 
intrinsic source T therefore vanishes in the bulk for this 
system. 

In the relaxation time approximation, the collision 
term in JSJ is given by 27\ : 



dt l 



fo 



where t p and t s are the momentum and spin relaxation 
times respectively, / is the identity matrix and a is the 
vector of Pauli spin matrices. In the extrinsic term in the 
source, the part depending on the momentum relaxation 
time will also cancel between the two bands, leaving us 
with just the contribution coming from the second term 
on the right hand side of (|13fl . The equation of continuity 
is then: 



dS_ 
~dt 



V-P r , 



(14) 



where P T = J d 3 kfp T is the torque dipole density. The 
two divergences will vanish in the bulk if the sample is 
homogeneous. 

We will now take a closer look at the spin current, 
making further use of the four-band model for the spin- 
orbit coupled valence bands of a zincblende semiconduc- 
tor. In previous work (ljj. Il4l llif . extensive discussions 
have been devoted to the extrinsic part of the spin cur- 
rent, which is given by the zero-field carrier velocity and 
spin integrated over the non-equilibrium part of the dis- 
tribution. Here we will concentrate on the intrinsic part 
of the spin current, coming from the field correction to 
the carrier spin current integrated over the equilibrium 
distribution. In order for this term to be dominant, scat- 
tering must be strong enough to keep the distribution 
function close to its equilibrium value, and weak enough 
to limit inter-band mixing. This is therefore opposite to 
the limit of Dyakonov-Perel |23| relaxation of spin in the 
weakly spin-orbit split bands of crystals. 

The two-fold degeneracy of both the heavy and light 
hole manifolds implies that, in the presence of an electric 
field, however weak it may be, mixing within the degen- 
erate manifold will occur in general. Fortunately, for the 
heavy holes the Jk = ±3/2 bands do not mix to first 
order in the electric field, and we can apply the nondc- 
generate band formalism. The s z dipole moment for the 
heavy holes is found to be = — ■ The torque 
dipole is, after an angular average, pjf = — e ^ z . The 
spin and torque dipoles are equal for both heavy hole 
bands. For the convective part of carrier spin current, in 
addition to a field correction to the carrier velocity due 
to the Berry phase jlflj . we obtain a term which is due 
to the change in the spin expectation value induced by 
the electric field and has the form iJ^E • ^j^-. Using 
these results, we find the current for spin-z component 
of a heavy-hole carrier to be 



3'h = ( 



1 



1 



1 



4fc 2 6m h A 12k 2 6k 2 



)eE x z, 



(15) 



•p 



where A = % — e; is the energy difference between the 
heavy and light holes. The first two terms come from the 
convective part due to field corrections to the carrier ve- 
locity and spin respectively. The third term comes from 
the rate of change of the spin dipole, while the last one 



4 



comes from the torque dipole. The heavy-hole carrier 
spin current can be simplified to: 



eE x 



1 



Jh — 



3fc 2 



(16) 



For the light holes, we must consider field induced mixing 
between the two degenerate bands. The details of this 
calculation will be deferred to a future publication [25| . 
we only quote the final result here which is very simple. 
The spin current per carrier in the light hole manifold 
has the same form as for the heavy holes, and differs 
only by a minus sign. Integrating over k and summing 
the contributions from all four bands, we arrive finally at 
the following expression for the total spin current: J s = 
csffE x z, where the spin-Hall conductivity is given by: 



e kh — h 



kh 



3tt : 



I _ JUL 

m h 



(17) 



where kh and ki are the Fermi wavevectors for the heavy 
and light holes respectively. Separate calculations based 
on the Kubo formula by the present authors and by Mu- 
rakami et al. |30j yield the same results. 

Finally, we comment on the relationship between bulk 
spin currents and spin accumulation near the edge of the 
sample. A theory of spin accumulation must start from 
the spin density continuity equation (|14fl . If the torque 
dipole density P T were absent from this equation, then in 
the steady state the spin accumulation would be due only 
to the spin current. The presence of the torque dipole 
density modifies the expression for the spin accumulation, 
giving that 



J Sdx = T s {J S x +Pl) 



(18) 



We have already discussed the response of the spin cur- 
rent to an electric field, and after a similar calculation 
for the torque density we find that 




mi 
2m,h 



-)• (19) 



Using n — 2.4 x 10 11 cm~ 2 and an electric field of 
20000V/cm as typical values, the spin current is — 10 25 
spins per unit area per second. We take the spin relax- 
ation time to be t s — 30ps 29] and a unit cell size of 
6.3A, and we obtain a spin accumulation of 1.2 x 10~ 4 
spins per unit cell area. This is a measurable effect as 
discussed in 0, 0] . 
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